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In one dimension, the exponential position operators introduced in a theory ol polarization are 
identified with the twisting operators appearing in the Lieb-Schultz-Mattis argument, and their 
finite-size expectation values zl measure the overlap between the unique ground state and an excited 
state. Insulators are characterized by Zoo 7^ 0. We identify zl with ground-state expectation values 
of vertex operators in the sine-Gordon model. This allows an accurate detection of quantum phase 
transitions in the universality classes of the Gaussian model. We apply this theory to the half-filled 
extended Hubbard model and obtain agreement with the level- crossing approach. 
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The metal-insulator transition is a fundamental phe- 
nomenon in strongly correlated electron systems. An 
insulator is distinguished from a conductor at zero tem- 
perature by its vanishing dc conductivity (Drude weight). 
Kohn argued that localization of the electric ground-state 
wave function is the signature of an insulating state [Q . 
Recently, Resta emphasized that insulators are charac- 
terized by their polarizability, and that meaningful defi- 
nitions are required for the polarization and position op- 
erators in extended systems. To this end, he discussed 
the ground-state expectation value of the exponential of 
position operators Xj in a finite size system |^-^, 
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He showed that the many-body expectation value of po- 
sition and polarization operators in periodic systems are 
related to ^Imlnz^. Generalization and numerical cal- 
culations of Zl were also made but our under- 
standing of this quantity is still far from satisfactory. 
Specifically, we need to know its relation to other the- 
oretical schemes, and to clarify the nature of phase tran- 
sitions detected by this quantity. 

In this Letter, we discuss the quantity zl from a differ- 
ent point of view. Limiting ourselves to one-dimensional 
(ID) cases, we give two interpretations to z^: One is 
based on the argument by Lieb, Schultz, and Mattis 
(LSM), which has been applied to spin and electron sys- 
tems to investigate structure of the excited states [0-0 . 
The other is the sine-Gordon theory p2|jr^ ], which de- 
scribes phase transitions in ID quantum systems in terms 
of the renormalization group. We will show that zl mea- 
sures the orthogonality between the unique ground state 
and a specific excited state of the finite size system, and 
that it also gives the ground-state expectation value of 
a vertex operator in the bosonization theory. Moreover, 
we show that the condition zl = corresponds to tran- 
sition points that belong to the universality class of the 



Gaussian model, and demonstrate this notion in numer- 
ical analysis for a lattice electron model. 

First, we discuss the new quantity based on a LSM- 
type argument using the ID spinless fermion model 
(SFM) 



^ = -tY.^4c^+l + H.c.) + V{{n,}), 



(2) 



where the number of lattice sites L is even, Ci is the 
fermion annihilation operator, and V{{ni}) denotes the 
interaction term given by arbitrary functions of the num- 
ber operator Ui = c|ci. The Fermi point is given by 
fcp = -kN/L where N is number of fermions. To make 
the ground state nondegenerate, we choose periodic (an- 
tiperiodic) boundary conditions for N = odd (iV = even), 
based on the Perron-Frobenious theorem. Applying to 
the ground state |\l/o) with momentum fc = the unitary 
"twisting operator" lA = exp[(27ri/L) ^^^j^ jn^] q times, 
generates a set of low- lying excited states \'^q) = U'^\'i>o). 
U is the lattice version of the operator appearing in Eq. 
(|l|). The translation operator T{cj — *■ Cj+i) has eigen- 
values e'*^. Due to the relation U-iTW ^ Te'^-f^F, the 
operator turns out to move q fermions from the left 
Fermi point (— fcp) to the right one {-\-k-p). The excitation 
energy for is evaluated as 

AE = {■^o\{U-mW - n)\^a) (3) 

L 

= -t[cos{2qn/L) - 1] Y.{^o\{cU+i + H.c.)|^'o>- 

Thus, if the state \'^q) is orthogonal to the ground state 
l^'o), there exists an excited state with energy oi 0{L^^). 

The orthogonality of these two states depends on the 
momentum of the excited state 2qkp. When fcp 7^ 7rp/g, 
where p is an integer and p/q is an irreducible fraction, 
1 5*0) and l^'q) are characterized by different quantum 
numbers, so that these two states are orthogonal. On 
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the other hand, when k-p = irp/q, the two states may not 
be orthogonal, so that this relation gives a necessary con- 
dition for the system to have a gap The overlap 
of the two states in a system of size L is given by 

z'f^ = (vl/olvf,) = {^„\W\^o)- (4) 

is real due to parity invariance (cj ~> c^+i-j). If the 
system is gapped (gapless), 2:00 7^ {zoo — 0) is expected. 
However, remains finite in finite-size systems even in 
gapless cases. The expectation value of the twisting oper- 
ator is nothing but the quantity underlying Resta's defi- 
nitions of electronic localization and polarization , and 
its generalizations ||^. 

Next, we consider zl in a different point of view. The 
low-energy excitations of the SFM (||) with kp = irp/q 
are described by the sine-Gordon model p^ , p^ 

nsG ^^J^ [Kid^ef + K-\d,cbf] 

+ 7^r^ [ da;cos[gV20(a:)], (5) 

with a relation [(t){x),d{x')] — — (i7r/2)sign(a; — x'). This 
effective model consists of the Gaussian model and the 
nonlinear term. Generally, the velocity u, the Gaussian 
coupling K, and the Umklapp scattering amplitude g 
are determined phenomenologically. The renormalization 
group (RG) equations for the nonlinear term are derived 
under a change of the cutoff a — > e^'a jl^ , 

where yo — 2{q^K/A — 1), j/i ~ q/t^v, and / = InL. 




FIG. 1. RG flow diagram for the sine-Gordon model. Phase 
transitions described by this diagram are the WZNW type 
(j/i =3/0 = 0), the BKT type = j/o > 0), and the Gaus- 
sian type {yi =0,yo < 0). The condition zl =Q corresponds 
to the case y\ = G which includes the WZNW-type and the 
Gaussian-type transitions. 



The RG flow of this equation shown in Fig. |^ de- 
scribes the following three type of transitions: As is well 
known, the SFM can be mapped on spin chains by a 
Jordan- Wigner transformation. When the system has 
SU(2) symmetry, Eq. (||) belongs to universality class 
of the SU(2)i symmetric Wess-Zumino-Novikov-Witten 
(WZNW) model |T^, and the parameters are fixed as 
2/1 = yo (9 = 2). Then, when the initial value of yi 
is positive yi(0) > 0, the nonlinear term is marginally 
irrelevant [j/i(cxd) = 0], and the system is gapless. On 
the other hand, for yi(0) < 0, the nonlinear term is rel- 
evant [2/1(00) —00], a gap opens, and the phase field 
is locked in the potential minimum as V8(0) = (mod 
2tt). We call the transition at yi(0) = "WZNW type". 
In the U(l) symmetric case where the ratio of ?/o and 
yi is no longer fixed, a Berezinskii-Kosterlitz-Thouless 
(BKT) -type transition takes place at |yi(0)| = yo(0) > 0. 
For the gapped states |yi(0)| > yo(0) > 0, the pa- 
rameters are renormalized as yi(oo) — ±cxd, and the 
phase field is locked as \/8{(j)) = tt, (mod 2tt). In 
the SFM at half-filling (spin chains with zero magnetic 
field), these two gapped states correspond to the charge- 
density-wave (Neel) and the bond-order-wave (dimer) 
states, respectively. On the unstable Gaussian fixed line 
[t/i(0) = with 2/0(0) < 0], a "Gaussian transition" takes 
place which is a second-order transition between the two 
gapped states, and the system is gapless on the transi- 
tion point. These three transitions can be identified by 
observing appropriate level crossings in excitation spec- 
tra in finite-size systems [^,^ . 

Now, let us interpret z^ in terms of the sine-Gordon 
theory. In the Tomonaga-Luttinger liquid theory which 
describes the Gaussian fixed point (j/i = 0), the cur- 
rent excitation created by the operator U corresponds to 
the vertex operator exp(i\/20) [|12| , p"3|jl7| . Besides, when 
kp — irp/q, the phase field changes as (/) — s- —4> under 
a parity transformation, so that zl given by Eq. (^) is 
related to the ground-state expectation values of the non- 
linear term as 

zi') = (cos(gV20)). (7) 

Since the sign of the nonlinear term in Eq. (^ changes 
at 2/1(0) = in the RG fiow diagram, the WZNW-type 
and the Gaussian-type transition points can be detected 
by observing z^ = 0. In the infinite-size limit, we expect 
Zoo = ±1 {zoc = 0) for the gapped (gapless) states. 

To demonstrate the above argument, we consider the 
ID extended Hubbard model (EHM) 

L 

Wehm = -tJ2Yl (4c»+i,. + H.c.) (8) 
+ U^^riiiTiii + niUi+i, 

i i 

at half-filling and zero magnetic field (kp = tt/2). Here 
L is even, and Cis is the electron annihilation operator 
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for spin s =t,i- The number operators are defined by 
= cl^Cis and rii = ni-^+riii. The effective Hamiltonian 
density of this system is given by sine-Gordon models for 
the charge {v = p) and the spin {v = a) sectors and a 
charge-spin coupUng term as 



v—p.cr 



(9) 



+ — TTT cos 



(27ra)2 
(27ra)2 



(27ra)2 
cos[V8(f>p] cos[\/8(/)(j] 



The spin part of this effective model belongs to the uni- 
versality class of the SU(2)i WZNW model. According to 
the level-crossing approach , the Gaussian transition 
in the charge part {g3± = 0), and the WZNW- type spin- 
gap transition {gi± = 0) take place independently near 
the U = 2V line with U > 0. Therefore, spin-density- 
wave (SDW), bond-charge-density- wave (BCDW), and 
charge-density-wave (CDW) phases appear that corre- 
spond to the locked phase fields {{(j)p) , {(pa)) = (0, *), 
(0,0), and (7r/\/8, 0), respectively, where * denotes the 
unlocked (gapless) case. 

To apply our argument to this electron system, we in- 
troduce the twisting operators for the charge and the spin 
sectors following Ref. pi] 



(10) 



where Ug = Gxp[{2ni/L)'^^^^ jujs]- Since this unitary 
operator corresponds to the vertex operator given by 
exp[i-\/2(0p + S(/)(j)], we obtain the expectation values of 
the nonlinear terms of Eq. (Q) as follows: 



=F(cos(\/8(/)i.)), V 



(cos(\/80p) cos(\/8(/)o-)). 



(11a) 
(lib) 



Therefore, (z^ 



P7 '^pa 



) are expected as (—1,0,0), 
(-1, 1, -1), and (1, 1, 1) for the SDW, BCDW, and CDW 
regions, respectively, in the L — > oo limit. We calculate 
these quantities in finite-size rings [L = 8-16) numer- 
ically. Eigenvalues and eigenstates of the EHM are ob- 
tained by the Lanczos algorithm and the inverse iteration 
method, respectively. 

In Fig. H, we show the numerical results near the 
U ~ 2V line at U/t = 3. For both charge and spin 
sectors, Zp and Zg- change continuously, and change their 
signs at different points: the Gaussian (BCDW-CDW) 
and the spin-gap (SDW-BCDW) transitions. Although 
the convergence of Zp, z^, and Zp^ to their saturation 
values is slow except for the transition points z^ — 0, 
they are expected to take the predicted values 0, ±1, 
and change discontinuously at the transition points in 
the L — > cx) limit. 
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FIG. 2. Behavior of Zp, Za-, and Zp^ defined in Eq. ( |l l| ) 
at U/t = 3 for L — 8, 10, 12, and 14 systems. In L ^ od 
limit, {zp, Zcr, Zpcr) are expected as ( — 1,0,0), ( — 1,1,-1), and 
(1, 1, 1) for the SDW, BCDW, and CDW regions. 
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FIG. 3. Size dependence of (a) the Gaussian transition 
point in the charge sector and (b) the spin-gap transition 
point at U/t — 3, obtained by monitoring z,^ — and the level 
crossings (I.e.). The system sizes are L = 8, 10, 12, 14, and 16. 
The extrapolation is done by K(i) = K(cx)) + A/L^ + B/i". 

In Fig. ||, we show the size dependence of the SDW- 
BCDW and the BCDW-CDW boundaries at U/t = 3 
determined by the conditions z^ — 0. A benchmark 
for their accuracy is provided by the level-crossing ap- 
proach which is also shown. These two methods give 
quantitatively similar results, so that our interpretation 
of the zl by the sine-Gordon theory is confirmed. For 
extrapolating the critical point to infinite system size, 
we use Vc{L) = Vc{oo) + A/V^ -\- B/L* given by confor- 
mal field theory, which is justified when the non-linear 
sine-Gordon term is absent. Numerical results equiva- 
lent to ours were also obtained in Ref. |l^ by observing 
discontinuities of Berry phases explained in the follow- 
ing. However, the physical interpretations are different. 
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Besides, the present calculation has great advantage in 
the accuracy and the computational time. 

We now discuss the quantity 2;^ as a complex variable. 
zl can be rewritten as a product of the discretized flux 
state 



M-l 



^m^{2qiT)) = n (vl/($,)|*(<i>,+i)), (12) 



J=0 



where $j = 2qjTT/M with M an integer. Here | *(<!>)) 
is the ground state for the Hamiltonian U^^{^)T-CU{^) 
withZ//($) = exp[(i$/L) ^^^-^ jrij]. Now we consider the 
phase angle of in the complex plane 7^ — Im Inz^. 
According to Resta ^ , this angle is regarded as the Berry 
phase |4||,|l|-|| 



2q7T 



d$(*($)|5ci.*($)). 



(13) 



The angle 7^ is or tt according to the signs of z^, so 
that it changes discontinuously at the Gaussian- or the 
WZNW-type transition point [yi(0) = with ?/o(0) < 0]. 
However, it also shows a discontinuity on the stable Gaus- 
sian fixed line [yi (0) = with yo (0) > 0] , which does not 
correspond to any transitions (see Fig.|^). In this case, jl 
becomes meaningless in the L 00 limit. For example, 
this situation appears in the charge sector of the EHM, 
near the U ^ 2V line with U < |l|jl|]. Extending 
Zl into the complex plane may be useful to distinguish 
phases with the same Re zl- For example, the SFM (^) at 
half-filling with V{{ni}) = +^'"■1^*4+2) have 

a charge-ordered state [{n2k) = (f^2/c-i-i) = [1 + (^l)''']/2] 
for V' ^ V > t that gives the same value z^' = — 1 
as that of the bond-order-wave state [{rij) = 1/2]. We 
speculate that their Berry phases 7^ are different, e.g. 
7i — — TT and jl = tt. It would be interesting to derive 
a simulation strategy (necessarily involving parity break- 
ing) to support such a conjecture. 

In summary, we have given the following two interpre- 
tations to Resta's expectation value of exponential posi- 
tion operators zl- One is as expectation values of a twist- 
ing operator which measures the orthogonality between 
the unique ground state and an excited state [Eq. (^]. 
The other is the ground-state expectation value of the 
nonlinear term of the sine-Gordon model [Eq. (^]. From 
the latter point of view, it is shown that insulating states 
characterized by Zqo = ±1 correspond to two different 
fixed points in the RG analysis, and that the condition 
Zl = gives phase transition point which belongs to the 
Gaussian universality class. We have demonstrated this 
notion in the EHM, and checked the validity of our ar- 
gument comparing with the results of the level-crossing 
approach. 

Our work could open significant extensions: The 
present argument can be applied to spin systems, with 
e.g. Haldane gaps and magnetization plateaus. The 



quantity zl may be applied to higher dimensional cases, 
since the LSM argument has been extended to these 
cases by considering a lattice wrapped on a torus |0,^,|5| . 
Besides, the quantity zl can be easily calculated nu- 
merically by the density matrix renormalization group 
method or quantum Monte Carlo simulations, as it is 
based only on ground-state properties. 
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(Heisenberg fellowship) and V0436/7-1. 
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